Strong spin-orbit coupling effects on the Fermi surface of Sr 2 Ru04 and Sr 2 Rh04 
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We present a first-principle study of spin-orbit coupling effects on the Fermi surface of Sr2Ru04 
and Sr2Rh04. For nearly degenerate bands, spin-orbit coupling leads to a dramatic change of the 
Fermi surface with respect to non-relativistic calculations; as evidenced by the comparison with 
experiments on S^RhCU, it cannot be disregarded. For S^RuCU, the Fermi surface modifications 
are more subtle but equally dramatic in the detail: spin-orbit coupling induces a strong momentum 
dependence, normal to the Ru02 planes, for both orbital and spin character of the low-energy 
electronic states. These findings have profound implications for the understanding of unconventional 
superconductivity in SnRuCU. 

PACS numbers: 71.18.+y, 74.70.Pq, 79.60.-i 



The emergence of unconventional superconductivity 
in Sr2Ru04 from a rather conventional Fcrmi-liquid-likc 
state has attracted the attention of numerous theoretical 
and experimental studies of this and related materials. 
Although Sr2RuC>4 is isostructural with the high-T c su- 
perconducting cuprates, there seems to be little doubt 
that the normal state is well described by a Fermi liquid, 
with a Fermi surface obtained from density-functional 
theory (DFT) band structure calculations 0, Q which 
agrees well with experimental determinations. The Fermi 
surface of Si~2Ru04 has been measured with several 
techniques with an unprecedented accuracy: de Haas- 
van Alphen (dHvA [|), Shubnikov-dc Haas (SdH d), 
low and high energy angle-resolved photocmission spec- 
troscopy (ARPES @, @|), and Compton scattering (CS 
0) data are all available for this compound and give a 
consistent outcome. The agreement of these experiments 
with calculations seems satisfactory but it is not a trivial 
result, if one recalls the large many-body renormalization 
of the band-dispersion Q and the notion that Sr2Ru04 
can be turned into a Mott insulator by Ca doping ||. 

The apparent much weaker influence of electron cor- 
relation effects in the Ad ruthenates, as compared to 
the 3d cuprates, is not unexpected. The radial extent 
of the 4c? wavefunctions with one radial node is much 
larger than the zero-node 3d one, leading to a larger one- 
electron bandwidth and a smaller on-site Coulomb re- 
pulsion. This fundamental difference between 3d and Ad 
transition metal compounds also explains the dominance 
of the t2 g band and low-spin-likc behavior in describing 
the low-energy scale physics for systems with less than 
6 d-electrons. In other words, Hund's rule splittings are 
eclipsed by crystal-field and band-structure effects in Ad 
systems. All this would suggest that Sr 2 Rh04, with one 
extra Ad-t^g electron, should also be well described by 
band theory and could provide an interesting approach 
to study the behavior of the very unconventional super- 
conducting state of Sr2Ru04 upon electron doping. Al- 



though one would expect the Fermi surface to change 
with the addition of one extra electron to the Fermi sea, 
it came as a surprise that band theory - within an ap- 
proach equivalent to that used for Sr2Ru04 - failed badly 
in describing the Fermi surface of Sr2Rh04, as deter- 
mined experimentally to a high degree of precision by 
both ARPES [13, EH and dHvA experiments 

The subject of this paper is to understand the rea- 
son for this discrepancy. It has been suggested that the 
failure of DFT in the case of Sr2Rh04 is a consequence 
of many-body interactions 13, 11 1 ■ This seems however 
an unlikely scenario since, on general grounds, many- 
body effects are expected to be of the same magnitude in 
Sr2Rh04 and Sr2Ru04. Indeed, as confirmed by experi- 
mental determinations of mass enhancement 0, E3| and 
quasiparticle dispersion [3, [13] , the renormalization pa- 
rameters are closely comparable for the two compounds 
and, at most, would suggest slightly larger renormaliza- 
tion effects for Sr2Ru04. We will show that electron 
correlations are not the main driving force here; instead, 
the problem can be well understood by taking into ac- 
count the interplay between spin-orbit coupling (SOC) 
and the details of the band dispersion close to the Fermi 
energy. In addition, the inclusion of SOC leads to a non- 
trivial, hitherto unknown, momentum dependence for or- 
bital and spin characteristics of the near-i?^ electronic 
states in Sr2Ru04; their knowledge is at the very heart 
of the microscopic understanding of spin-triplet, possibly 
orbital-dependent, superconductivity [l3j |. 

Experimentally it was shown that SOC does play an 
important role in insulating ruthenates 14|. This is not 



surprising since the atomic relativistic SOC constant is 
C=161 (191) mcV for Ru 4+ (Rh 4 +) [IIEl, which is not 
a small quantity. Furthermore, it is also rather similar for 
the two systems leading one to believe that, if SOC is not 
important for the metallic ruthenates, the same should 
hold for the rhodates and viceversa. As we will show, the 
inclusion of SOC dramatically improves the agreement 



2 



Sr 2 RhQ 4 



Sr 2 Ru0 4 




Z-dispersion 



*4 



x4 



M 



-M 



M 



M 



FIG. 1: (color online). (a,b) LDA Fermi surface of S^RhCU and (c,d) Sr2Ru04; (e-h) k z electronic dispersion for the Sr2Ru04 
cuts highlighted in (c,d) by green solid bars. Calculations were performed without SOC ( top panels: a,c,e,g) and with SOC 
(bottom panels: b,d,f,h). The grey-scale ARPES data are reproduced from (a,b) Kim et al. [Tit and (c,d) Damascelli et al. [jj. 



between DFT and experimental results for S^RhC^; we 
will also explain why SOC has a much smaller appar- 
ent effect for Sr2Ru04. Care should be taken, however, 
since for Sr 2 Ru04 there are specific regions in the three- 
dimensional Brillouin zone where SOC does have a very 
large effect; this might have an important influence on 
the interpretation of several experimental findings, from 
the details of the k z electronic dispersion 0, H| to the sig- 
natures of spin-triplet superconducting pairing [I?], EH ■ 

DFT calculations were done with the code WIEN2k 
[lit and structural data reported in Ref.[2Cj|. To test the 
generality of the results, we used several values for the 
radius of the muffin-tin spheres (RMT) and two differ- 
ent functionals, local density approximation (LDA) and 
generalized gradient approximation (GGA). The results 
are basically independent of RMT and extremely simi- 
lar for LDA and GGA, especially close to the Fermi en- 
ergy. The plots and values shown in the paper were calcu- 
lated within LDA with RMT Sr =1.9a , RMT flu =1.75a , 
RMTo=1.45ao. We used a cutoff parameter RK max = 
9.0; eigenstates with energies up to lORyd were included. 
Self-consistency was obtained with ~400 and ~100 k- 
points in the reduced Brillouin zone for Sr 2 Ru04 and 
Sr2Rh0 4 , respectively. The Fermi surface was calculated 
on a square fc-mesh, with points separated by -^tt within 
the plane of the calculated surface. 

In Fig.QJ, and[TJ; we show the Fermi surface of Sr2Rh0 4 
and Sr2Ru04 calculated without SOC, which is in good 
agreement with previous theoretical results 0, 0, EH 11 1 . 
The Fermi surface of Sr2Ru04 consists of three sheets, 
labelled a, (3 and 7. The 7 sheet is mainly composed 
of d xy orbitals and is highly two dimensional; the a and 
(3 sheets are formed by the one-dimensional-like bands 
derived from d xz and d yz orbitals, which exhibit anti- 
crossing behavior (i.e. mixing) along the zone diagonal. 
In Fig.[TJ; we present calculations for both k z = and 
7T (centered around T and Z, respectively), which show 



some degree of k z dispersion. As k z is a difficult to deter- 
mine quantity in the experiment the comparison be- 
tween LDA and ARPES should be attempted for various 
values of k z . For k z = 7r, the agreement between ARPES 
H, dHvA Q (not shown), and LDA is rather good; for 
k z = 0, however, the crossing between 7 and (3 sheets in 
LDA is experimentally not reproduced. Note that this is 
really a crossing and not an anticrossing, as for k z = 
the bands arising from d xy and d xz /d yz orbitals are of 
different symmetry and no mixing is allowed. Also, con- 
trary to the LDA results of Fig.[T}:, quantum oscillation 
experiments 0, Q] indicate very little k z modulation for 
the Fermi surface (smaller than 1% of the zone even for 
the (3 sheet, which exhibits the largest k z effects). 

In principle the same discussion holds for Sr 2 Rh04; 
there are however a few complications due to structural 
distortions. The RhOg octahedra in Sr 2 Rh04 are rotated 
around the c-axis, which results in a y/2 x y/2 doubling of 
the unit cell in the a-b plane, and also a doubling along 
the c-axis because of the alternation in the RhOe rotation 
direction. This leads to a doubling of all bands and Fermi 
sheets, as evidenced in Fig.QJt by the back-folding of the 
Fermi surface sheets with respect to the line connecting 
two closest M points, e.g. (7r,0,0) and (0,7r,0). In ad- 
dition, this distortion reduces the electronic bandwidth 
and allows mixing between d xy and d x i_ y i orbitals ll|, 
which causes the 7 sheet to almost disappear; small re- 
maining 7 pockets can be seen around the a and (3 sheet 
crossings. The agreement between LDA and experiment 
for the rhodate is, as shown in Fig.QJi and also noted in 
the literature 13, 111, not very good (especially if com- 



pared with the beautiful agreement found for Sr 2 Ru04). 

At present, there is no explanation of why LDA is able 
to describe Sr2Ru04 but fails for Sr2Rh04, two materi- 
als that are structurally and electronically quite similar. 
We will show here, starting from the more striking case 
of Sr2Rh04, that a resolution of this discrepancy can be 



obtained with the inclusion of SOC in LDA. As one can 
see in Fig.QJi, the agreement between LDA+SOC and 
experiment is very satisfactory for S^RhC^: the 7 sheet 
is now completely absent, and a and (3 sheets are sig- 
nificantly smaller (the a pocket is still somewhat larger 
than in the experiment, which could possibly be solved by 
fine-tuning the Ru06 distortion at the surface [ll[ but is 
beyond the point here). As for Sr 2 Ru0 4 (Fig.EEH), SOC 
effects are much smaller, especially for k z = 7r; some im- 
provement can be observed only for k z = 0, in which case 
the 7 and (3 sheet crossing becomes an anticrossing be- 
cause the bands no longer have different symmetry. This 
point is important, as we will stress below, but first we 
discuss the reason why the effect of SOC for Sr2Rh04 
seems so much stronger than for Sr2Ru04. 

Fig-El presents the band structure of Sr 2 Ru04 and 
Sr 2 Rh0 4 from LDA (top) and LDA+SOC (bottom). 
It is clear that if the i 2g bands are already well sep- 
arated in energy, e.g. at the M point for Sr 2 Ru04 in 
Fig. UK SOC only mildly shifts the LDA eigenstates. A 
whole new energy splitting is instead induced by SOC at 
those momenta characterized, in LDA, by degenerate £ 2g 
bands. For instance in Fig. [2}:, the degenerate dxi^ and 
dyl"^ LDA bands of Sr 2 Ru04 along T-Z are recombined 
into two doubly-degenerate complex bands ,c?L 1 and 
d l _ 1 ,d\ with d±i = ^^(d xz ± idyz)], split by an energy 
equal to the SOC constant £ = 93 meV. Note that the 
latter is derived from the splitting at the T point, and 
the strong reduction with respect to the 161 meV atomic 
value of Ru is mainly due to the dilution effect resulting 
from the covalent mixing with the 0-2p bands (accord- 
ingly, one also finds a SOC splitting for the 0-2p derived 
bands around -3.8 eV due to the mixed-in Ru character; 
and the SOC splittings are larger at X and M points, 
where the covalent mixing is smaller). One might won- 
der, in this regard, whether the larger SOC effects in 
Sr 2 Rh04 might arise because of the ~ 20% larger atomic 
SOC for Rh (191 meV) as compared to Ru, originating 
from the somewhat less extended Ad orbitals. However, 
the on-site energy difference for Rh-4c? and 0-2p orbitals 
is smaller than for Ru-4d and 0-2p, making Sr 2 Rh04 
more covalent than Sr 2 Ru04. As a result, the effective 
SOC for the i 2g bands in the two compounds is about 
equal (Fig.[2j:,d); the magnitude of the atomic SOC con- 
stants is thus not the reason for the different behavior. 

More clues on the magnitude of SOC effects for the two 
compounds come from the inspection of the Sr 2 Rh04 
band dispersion in Fig.[2j},d. The structural distortion 
responsible for band folding and narrowing in Sr 2 Rh04 
gives rise to several nearly-degenerate states close to the 
Fermi level (e.g., the d xy band maximum at X in Fig.[2]a, 
is back- folded to the now equivalent L point, leading to 
the mixing with the e g band and the appearance of two 
da;j,-derived maxima just above Ep, between T and X 
in Fig. (2b). Because of the presence of weakly-dispersive 
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FIG. 2: (color online). (a,c,e) Sr 2 Ru0 4 and (b,d,f) Sr 2 Rh0 4 
LDA dispersion calculated without SOC (a,b) and with SOC 
(c-f). (e,f) Enlarged view of the marked i^F-regions in (c,d). 



bands close to Ep, the SOC-induced splitting causes sub- 
stantial modifications of the Fermi surface in Sr 2 Rh04 
(yet in the full respect of Luttinger's theorem), such as 
the complete disappearance of the 7 sheet (Fig.[2f). On 
the contrary in the case of Sr 2 Ru04 in the very same k- 
space region (Fig. [2k), the energy shift from SOC does not 
translate in large changes in k because of the steep band 
dispersion. It is therefore the interplay between SOC and 
band filling (i.e., the presence, or lack thereof, of weakly- 
dispersive degenerate states close to the chemical poten- 
tial), which determines the more dramatic Fermi surface 
changes for Sr 2 Rh04 in LDA+SOC. Larger effects should 
be expected for Sr 2 Ru04 upon increasing the electron fill- 
ing; this has indeed been observed in Sr 2 _ y Laj,Ru04 [2l| . 
in which the 7-sheet topological transition detected for 
y>0.20 is well accounted for by LDA+SOC (not shown). 

In the following, we will focus on some aspects of the 
LDA+SOC results which, beyond the detailed analysis 
of the Fermi surface, promise to have profound implica- 
tions for the physical properties of Sr 2 Ru04. We have 
already noted, in relation to Fig.[TJ:, that for k z = the 
LDA calculations display the crossing of 7 and (3 sheets 
at eight momenta symmetrically located with respect to 
the T-X lines. As SOC is not a small perturbation at 
highly-degenerate /c-points in LDA, to understand the 
impact of SOC on the low-energy physics we have to fol- 
low this degeneracy as a function of k z . In Fig.[T^,g we 
plot the k z modulation of the LDA Fermi sheets in the 
T-X-Z and T-M-Z planes. The fc-space loci of 7 and 
(3 sheet degeneracy in LDA defines an elliptical contour 
centered around the T-X line, crossing the k z = plane 
at (k x ,k y )=(0.28Tr,0.34ir) and (0.347r,0.287r) as indicated 
in Fig. [Hi, and crossing the T-X-Z plane at k z ~ O.6871". 
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This is a degeneracy between the d xy orbital and a lin- 
ear combination of d xz and d yz orbitals; there is also a 
further degeneracy at (k x ,k y ,k z ) = (0.3Air,0.3Air,0.33ir) be- 
tween the d xz and d yz orbital (Fig. []};). In addition to lift- 
ing the degeneracy of the LDA Fermi surfaces (Fig.[Tf), 
the inclusion of SOC reduces to more realistic values the 
k z - modulation of each individual Fermi surface sheet, 
which in LDA is much larger than observed in quan- 
tum oscillation experiments 0, 0|. Note, in this regard, 
that the rather localized bulging of the Fermi surfaces 
in LDA+SOC, such as e.g. along the zone diagonal for 
k z = for the 7 sheet, might require the inclusion of 
higher order cylindrical harmonics than usually consid- 
ered in the quantum oscillation analysis 0, 0] • 

Perhaps more important, SOC is responsible for the 
three-dimensional fc-dependence of the wavefunction spin 
and orbital character, along the Fermi surface. For in- 
stance, the character of the 7 sheet is no longer pure d xy ; 
instead, for k x = k y and \k z \ < 0.687T, it is mainly a 
linear combination of d xz and d yz . This not only influ- 
ences the hopping parameters; it also has a large effect 
on the spin anisotropy. At (0.347r,0.347r,0.337r), SOC re- 
combines the degenerate d xz ^ and d$ z '^ orbitals into 
the two doubly-degenerate gJx,gJ1i and d^_ 1} d\ complex 
orbitals. These new eigenstates have a spin- moment in 
the z direction and a zero expectation value for x and 
y spin components. This follows from s x = \{s + + s~) 

and (djls - ^^) = 0, as the orbitals are orthogonal to 
each other, and analogous relations for s + and s y . Thus, 
at (0.347r,0.347r,0.337r) the spins on the 7 and a sheets 
point along the c-axis, as indicated in Fig.[TJi. Similar 
arguments can be made if the d xy and d xz or d yz or- 
bitals are degenerate. By cyclic permutation of coordi- 
nates, it can be shown that this leads to a spin direc- 
tion in the a-b plane (Fig.Q}:), along the whole ellipti- 
cal LDA-degeneracy contour between 7 and (3 sheets dis- 
cussed above. Note that this fc-dependent non-collincar 
spin anisotropy cannot easily be seen in the bulk mag- 
netic susceptibility, which measures an average over all 
fc-points. It has however severe implications for all fc- 
dependent properties of Sr2Ru04. 

Our study shows that, when describing experiments on 
Sr2Ru04 and Si'2Rh04, it is crucial to start from accu- 
rate ab-initio band-structure calculations including SOC. 
One should be careful about a phenomcnological mod- 
elling of the electronic structure. For instance, a simple 
tight-binding fit to the experimentally determined Fermi 
surface of Sr2Ru04 leads to non-degenerate bands at the 
Fermi energy; adding SOC, starting from such a scenario, 
would only introduce a small perturbation. SOC is how- 
ever not a small perturbation for metals if (nearly) de- 
generate bands cross or are close to the Fermi energy, 
which happens for Sr2Rh04 in extended portions of the 
Brillouin zone and for Sr2Ru04 at some momenta. This 
has important implications - macroscopic in the case of 



Sr2Rh04 - for fundamental properties such as the de- 
tailed shape and orbital character of the Fermi surface 
sheets, as well as the amount of k z dispersion. For exam- 
ple, SOC is responsible for a significant reduction of the 
effective three-dimensionality in Sr2Ru04, which is con- 
sistent with the dHvA results Also, even in the over- 
all less dramatic case of Si - 2Ru04, SOC leads to a spin 
anisotropy and non-collinear behavior in several fc-space 
regions, with the spin precession direction varying with 
k according to the orbital character of the LDA+SOC 
eigenstates. This fc-dependent orientation of the expec- 
tation value of the spin calculated on the LDA+SOC 
eigenstates, which reminds us of Rashba coupling, could 
be a very important ingredient in describing the uncon- 
ventional superconducting state of Sr 2 Ru0 4 and should 
be investigated further. Since all Fermi surface degenera- 
cies are lifted in LDA+SOC and the fc-dependent energy 
splittings are much larger then the superconducting pair- 
ing energy, the SOC-induced spin anisotropy in concert 
with orbital mixing should directly affect the orbital and 
spin angular momentum of the Cooper pairs. In fact, 
singlet and triplet states could be strongly mixed, blur- 
ring the distinction between spin-singlet and spin-triplet 
pairing [l3| , and with it the simple interpretation of time- 
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reversal symmetry-breaking experiments 
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